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Using the criteria of simplicity, positivity, smoothness and process-independence of construction, 
we find a non-polynomial model for the quark distribution amplitude of the nucleon that is in agree- 
ment with QCD sum-rule moments. This simpler ("haplousterotic") model solves some of the prob- 
lems connected with the longitudinal momentum dependence of the wave function, whereas the k±- 
dependence of the wave function, which is also important, is still left undetermined. The expansion 
properties of the model into a polynomial series are examined. 



The quark wave function of the nucleon is of out- 
standing theoretical importance because it is a universal, 
process-independent quantity. However, the wave func- 
tion cannot be obtained directly by solving the equations 
of motion because of the complicated structure of QCD. 
A simplification of the problem can be made by treat- 
ing the transverse momentum dependence of the wave 
function in an effective way so that only longitudinal mo- 
menta have to be taken into account explicitly. This is 
done e.g. in the definition of the so called nucleon quark 
distribution amplitude <f>^ via the light-cone wave func- 
tion Ipx: 

<f> N (xi, n 2 ) := / [d 2 k u ] ip N (xi, k±i) . (1) 

<I>n(xi ; x 2, £3; A* 2 ) can be interpreted as the probability 
amplitude for the hadron to consist of quarks collinear 
up to the scale fi 2 with longitudinal momentum fractions 
Xi. tpN(x,k±) is the so-called light-cone Fock-state wave 
function of the hadron, which has to be understood as 
the probability amplitude for the hadron to consist of 
3 quarks with longitudinal momentum fractions and 
transverse momenta k 

Some models for the quark distribution amplitude <f>^ 
were used to construct the full wave function by mul- 
tiplying it with an additional function that has also a 
dependence on the transverse momenta, e.g. 

ip N (x, k L ) = <f> N {x) ■ f(x, k L ) . (2) 

Although the problem of constructing a realistic wave 
function for the nucleon in a reliable way is of outstand- 
ing interest, this task is far from being completed. On 
the contrary, a large number of different models for the 
nucleonic quark distribution amplitude has been put for- 




FIG. 1. Some representative models for the nucleon wave 
function appearing in the literature. The asymptotic form 
4>N, as = I2OX1X2X3 is given for comparison. The different 
amplitudes were taken from the following references: </>coz: 
[3]; 4>cz: [2]; 4>ks- [5]; 4> GS : [4]; het .: [8]; 4>™ , 4> 1Vc : [6]; <t> { $: 
[7]. All these models are polynomial models of 2 nd and 3 rd 
order. For a discussion on polynomial models see [9]. 



ward in the past by various authors. As the most im- 
portant examples, we show in Fig. 1 some models for the 
nucleon from Refs. [2-8]. Among these, the most recent 
one is the "heterotic" quark distribution amplitude [8] 
which was constructed to amalgamate characteristic fea- 
tures of the preceeding models of Refs. [4,3]. We see large 
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variations among the different wave functions which indi- 
cates that moments calculated from QCD sum rules are 
not sufficient in order to fix the shape of the distribution 
function in a reliable way. 

Recently, it has been shown in Refs. [9,10] that it is 
not possible to find polynomial models for hadronic dis- 
tribution amplitudes from their corresponding moments 
calculated by means of QCD sum-rules. The reason for 
this is the extreme instability of the polynomial expan- 
sion with respect to small variations in the moments. 
Already small deviations from the exact moment values 
result in a "neurotic" polynomial expansion, i.e. strong 
oscillations occur that completely distort the shape of the 
original function. This is discussed in detail in Ref. [10]. 
The experience with polynomial models indicates that 
oscillating wave functions are no serious candidates for 
the true distribution amplitude, even if additional con- 
straints are imposed [6,8]. 

However, for perturbative QCD calculations of exclu- 
sive processes it is necessary to have at least some simple 
model for the distribution amplitude to deal with. 

Consequently, one has to formulate some strong crite- 
ria that a candidate for a realistic model wave function 
has to fulfill in addition to QCD sum-rule moments so 
as to rule out the pathological cases [1,11-13]. An espe- 
cially important criterion is the requirement of simplic- 
ity. This means that the functional form of the model 
has to be as simple as possible with a minimum number 
of free parameters that can be adjusted to the moments. 
Furthermore, the sum-rule moments of Ref. [3] are pos- 
itive or allow positive values within their error bars. It 
is therefore sensible to require the function to be strictly 
positive. Another point is that, since polynomial expan- 
sions have turned out to be more than dangerous, the 
wave function has to differ essentially from a finite poly- 
nomial. Finally, the wave function has to be determined 
in a process-independent way, in contrast to the models 
in the past which were adapted to experimental values of 
form factors or decay widths. 

We summarize the criteria for the model distribution 
of the nucleon that are essential in our opinion: 
(i) functional simplicity, (ii) minimum number of free 
parameters, (iii) smoothness, (iv) absence of oscillations, 
(v) strict positivity, (vi) substantially non-polynomial 
(e.g. exponential) form and (vii) process-independent 
construction. 



The haplousterotic model 

According to these criteria, we look for a model wave 
function of simpler ("haplousterotic") type by modelling 
the completely symmetric asymptotic distribution am- 
plitude, mainly by shifting the position of the maximum. 
We want to stress that, if this leads to a function that 



fulfills all requirements, there would be a convincing al- 
ternative to other wave functions that have a more com- 
plicated structure. Since we know that the method of 
constructing a finite polynomial expansion directly from 
sum-rule moments is inapplicable, a possible candidate 
for such a wave function is an exponential function that 
is also smooth and positive by construction. We mention 
that functions of such type have already been suggested 
by the authors of Ref. [14]. 

As a specific ansatz for our "haplousterotic" (Ha + ) 
model of the quark distribution amplitude of the nucleon, 
we choose the following form: 

( b (r) b (r) b (r) \ 
4>*- + (x) = Nexp "(^ + ^ + ^) • (3) 

\ - L \ - L 2 ^3 J 

(r) 

The three parameters & s - are adjusted in order to change 
the position of the maximum according to the require- 
ments of QCD sum-rule moments, r = (-|, 1, 2, . . .) for 
model (A, B, C, . . .), and N is the normalization fac- 

(r) 

tor. In other words: the b { are chosen so as to yield 
optimal agreement between the moments of the model 
distribution amplitude, 

(ijk) := [x\x^ 2 x\) := J [dx] x\ x? 2 x\ ■ (f)f 1 a+ (xi, X2, £3) , 

(4) 

and the corresponding values of the same moments as 
obtained from QCD sum rules. 

We verify that our criteria for the wave function are 
indeed satisfied: (i) functional simplicity: exponential 
function, (ii) minimum number of parameters: the b^p 
are three independent parameters that are indispensable 
to determine the location of the maximum and its width, 
the model (iii) is smooth, (iv) has no oscillations, (v) is 
positive, (vi) substantially non-polynomial and (vii) no 
experimental information has been used to construct the 
model — the only input is given by QCD sum-rule mo- 
ments! 

For r = 1/2 model A reproduces the 19 moments of 
Ref. [3] within an accuracy of about 20 % with only three 
parameters! The corresponding moments are shown in 
Table I in comparison with those of <f> as Ref. [5] and 
Ref. [3]. Furthermore, also functions for different r show 
a strong similarity with model A, if the moments are suf- 
ficiently reproduced. We have therefore another criterion 
for the credibility of a model constructed from QCD sum- 
rule moments: the shape of the resulting model function 
should not depend too strongly on the specific ansatz 
used [12]. 

As an example, we show plots of model A in Fig. 2. 
In Fig. 3 we illustrate how the position of the maximum 
is shifted from its (asymptotic) position in the middle of 
the triangle. 
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TABLE I. Two different sets of sum-ruie moments of King and Sachrajda [5] and of Chernyak, Ogloblin and Zhitnitsky [3], 
in comparison with the moments of the following wave functions: </>jv,as, </>coz [3] and our new model (j> 1 ^ a+ . 

Moments (n\ «2 nz) 



ni n 2 n 3 




KS 


COZ 




(j>N,a,s 


0COZ 









1 


1 




1 


1 


i 


1 


0.46 


. . . 0.59 


0.560 


+ 0.06 
-0.02 


0.333 


0.579 


0.534 


1 
1 


0.18 
0.22 


. . . 0.21 
. . . 0.26 


0.192 
0.229 


+ 0.008 
-0.012 
+ 0.021 
-0.029 


0.333 
0.333 


0.192 
0.229 


0.209 
0.257 


2 


0.27 


. . . 0.37 


0.350 


+ 0.07 
— 0.03 


0.142857 


0.369 


0.323 


2 


0.08 


. . . 0.09 


0.084 


+ 0.004 
-0.019 


0.142857 


0.0680 


0.0704 


2 


0.10 


. . . 0.12 


0.109 


+ 0.011 
-0.019 


0.142857 


0.0890 


0.0956 


1 1 


0.08 


. . . 0.10 


0.090 


+ 0.01 
-0.01 


0.095238 


0.0970 


0.0942 


1 1 
1 1 


0.09 


. . . 0.11 


0.102 
-0.03 ... 


+ 0.008 
-0.012 

0.03 


0.095238 
0.095238 


0.113 
0.0270 


0.117 
0.0441 


3 
3 

3 
2 1 
2 1 

1 2 
1 2 
2 1 






0.236 
0.032 
0.052 
0.045 
0.050 
0.035 
0.041 

-0.004 ... 


+ 0.014 
-0.026 
+ 0.008 
-0.004 
+ 0.004 
-0.004 
+ 0.004 
-0.004 
+ 0.005 

— 0.006 
+ 0.002 

— 0.008 
+ 0.002 
-0.004 

0.007 


0.0714285 

0.0714285 

0.0714285 

0.035714285 

0.035714285 

0.035714285 

0.035714285 

0.035714285 


0.245 

0.0381 

0.0485 

0.0587 

0.0658 

0.0243 

0.0331 

0.0056 


0.212 

0.0300 

0.0435 

0.0489 

0.0617 

0.0277 

0.0382 

0.0127 


1 2 






-0.005 ... 


0.008 


0.035714285 


0.0073 


0.0138 


1 1 1 










0.0238095 


0.0141 


0.0176 



-120° 0° 120° 

* 




1 2. The typical shape of the "_ffa + " (simpler) distri- 
bution amplitude determined to reproduce the QCD sum-rule 
moments of COZ. agrees also with the QCD sum-rule 

moments of KS. Note that this model function has only three 
free parameters, which were determined only by a best fit to 
the COZ moments. No attempt was made to adjust form 
factors or decay widths. Thus, the parameters were deter- 
mined process-independently. The results do not depend very 
much on the precise form of the function (e.g. r = i,l,2). 
Although the function seems to enter the x 2 = 1-side very 
steeply, it reaches x t = with slope - as is clearly seen on 
the x\ = 0-side. For reference also V123 = §-(</>2i3 + ^123) and 
A 12 s = |(02i3 - ^123) are given. 




x 2 = l x 3 = l 



FIG. 3. This figure illustrates the deformation of the distri- 
bution amplitude away from its asymptotic form. Essentially, 
the maximum of the distribution amplitude moves from the 
completely symmetric starting point (®) very near to the side 
with X2 = and towards x\ = 1 (ffi). We found for several 
models always xf ax > 2/3. This simple method suffices to 
reproduce the QCD sum-rule moments within a reasonable 
accuracy. 



3 



It seems therefore that some of the problems of con- 
structing a (purely longitudinal) quark distribution am- 
plitude of the nucleon have been solved: the longitudinal 
momentum dependence of the wave function can now be 
fixed sufficiently by sum-rule moments without having 
recourse to experimental input. It follows that the main 
uncertainty in the construction of the full wave function 
lies in the still undetermined fcj_-dependence. 



Evolution and convergence properties of the model 

Since we have now an exact functional form for our 
model distribution amplitude of the nucleon, there is in 
principle no reason to use a polynomial expansion in 
practical calculations. However, in order to obtain the 
correct Q 2 -evolution of the distribution amplitude, the 
function has to be expanded into an infinite series of Ap- 
pell polynomials Ai{x) according to 

= ^ s (x) £ c * f^KvT' m*) ■ (5) 

i = \ a S(i l o) / 

As the precise functional form of our model is known by 
definition, see Eq. (3), the expansion coefficients c; can 
be determined directly via the orthonormality relations 
of the Appell polynomials: 

Ci = J [dx] Ai{x) <f,% a+ (x,n 2 ) . (6) 

We note that it is not necessary to use the Appell poly- 
nomials as constructed by Lepage and Brodsky [15] if the 
evolution effect is ignored. <j)fj a (x,/j'q) can then be ex- 
panded into the polynomials originally found by Appell 
[16]. Since the expansion is unique, it does not depend 
on the set of basis polynomials. 
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FIG. 4. The polynomial expansions of the Ha distribu- 
tion amplitude (cf. Fig. 2) in different orders. The coefficients 
of the Appell expansion of the best fit of </>jJ a+ are given in 
Table II. A polynomial expansion is needed if the evolution 
of the wave function is taken into account non-numerically. 
Although it seems that already 10* order resembles the wave 
function sufficiently, this is not the case for all functions of 
"iia + " type. 



It is demonstrated in Fig. 4 that the convergence of 
the whole series is quite slow: a very large number of 
polynomials has to be summed up until the shape of the 
exponential function is sufficiently reproduced. We see 
that, although the first order gives already a qualitative 
impression of the rough structure of the function, the ap- 
proximation is completely insufficient: the characteristic 
location of the maximum has moved far inside the tri- 
angle. We see only a very gradual improvement of this 
fact in the next orders of expansion, whereas seizable os- 
cillations appear (we recall that we are dealing with the 
exact polynomial expansion). In seventh and higher or- 
der the expansion becomes smoother and the maximum 
of the distribution amplitude is located in the end-point 
region as required by the exact function. It seems that 
an expansion of order ten or more is necessary to yield an 
approximation of the exact function that can be used in 
numerical calculations. The corresponding expansion co- 
efficients up to fifth order with respect to the normalized 
Appell polynomials of Lepage-Brodsky type are shown in 
Table II. We see that the coefficients remain of the same 
order of magnitude even for order 5. The slow conver- 
gence can thus be seen also numerically. 

TABLE II. Expansion coefficients of the model distribution 
amplitude c/>^ a up to fifth order with respect to the Appell 
polynomials of Brodsky and Lepage [15]. The polynomials up 
tp third order are arranged according to Ref. [7] and the re- 
maining ones according to Ref. [18]. The Appell polynomials 
were taken normalized to 1. 



order 


n 










1 


1 


1 


0.82916 




2 


-0.68548 


2 


3 


0.42216 




4 


-0.34454 




5 


0.26269 


3 


6 


0.19964 




7 


0.19908 




8 


0.17937 




9 


0.33544 


4 


10 


-0.20640 




11 


-0.12029 




12 


0.16415 




13 


0.25620 




14 


-0.04685 


5 


15 


0.19552 




16 


-0.02693 




17 


-0.13538 




18 


-0.05044 




19 


-0.17528 




20 


0.01135 



Summarizing our results we recall that using the crite- 
ria of simplicity, smoothness, positivity etc. one is able 
to construct a quark distribution amplitude which is in 
satisfactory aggreement with QCD sum-rule moments. 
We emphasize that, in view of the fcj_-dependence of the 
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full wave function, we did not adapt the distribution am- 
plitude to form factors or decay widths because the 
dependence leaves still a large freedom in the calculation 
of physical observables: In the form factor calculations 
the end-point regions give large contributions as has been 
shown by Isgur and Llewellyn-Smith [17]. Therefore, the 
^-dependence of the full wave function will have a sig- 
nificant effect on the form factor. The effects of this 
dependence have been studied by means of a mass in the 
gluon propagator or by taking into account Sudakov form 
factors [19,20]. Since the resulting values are sensitive to 
the specific choice of the cutoff or cutoff-procedure, there 
is still a large range of possible values for the form factor. 

Because of these reasons, we tried to determine the 
quark distribution in a process-independent way. To this 
end, we chose an ansatz for the model that is as simple as 
possible and adapted it to QCD sum-rule moments only. 
This model for the longitudinal momentum distibution 
can then be kept essentially fixed and further improve- 
ments can be attempted mainly in the construction of 
the correct fcj_-dependence of the full wave function. 

Finally, we would like to point out that because of 
its process-independence the quark wave function of the 
nucleon is an interesting object by itself, apart from per- 
turbative QCD calculations. The calculation of form fac- 
tors and other observables has nothing to do with con- 
structing the wave function. We believe that our "hap- 
lousterotic" model is the first simple model of the nu- 
cleon distribution amplitude that is consistent with the 
requirements of QCD sum-rule moments and for which 
convergence is explicitly verified. 
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